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ABSTRACT : 



An analysis for the transient axial response of a gun launched motor 
case with unbonded propellant is performed. The present effort extends a 
previous work by including higher harmonics in the representation of the 
breech force at the aft end of the motor case. Including the higher har- 
monics leads to a more accurate representation of the actual breech force 
The results of the analysis indicate that the dynamic effects on the in- 
ternal axial force in the case are negligible for a breech force duration 
of 0.026 seconds. This is the same conclusion that was made in the pre- 
vious work based upon a simple representation of the breech force. 
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area of motor case cross section 
Fourier cosine coefficients of breech force 
Fourier sine coefficients of breech force 
speed of wave propagation 
axial strain 

Young's modulus of elasticity 
axial force 
constant 
1 ength 

mass of motor case 
breech force 
time 

deformation 

displacement 

rigid body displacement 

transformation variable 

transformation variable 

coordinate axes 

coordinate axes 

time duration of breech force 

constant 

mass density of the motor case material 

axial stress 

constant 



INTRODUCTION 



This report examines the transient axial response of a gun-launched 
rocket motor, with unbonded propellant, during the launch phase. It is 
part of a continuing effort to develop structural design procedures for 
gun-launched rockets. In a previous report (Ref. 1), an analysis for the 
time-dependent axial force distribution in the motor case was performed 
in which the force at the breech end of the case was approximated by a 
one-half sine wave over the duration of the loading. That analysis estab- 
lished that for the one-half sine wave input, a time duration of .012 
seconds leads to negligible dynamic effects, i.e., for a one-half sine 
wave of .012 seconds duration, the internal axial force at any location 
along the case appears as a duplicate of the breech force, but of smaller 
magnitude. It was shown, however, that decreasing the time of the force 
duration to .0012 seconds leads to significant dynamic effects. This re- 
sult indicated that higher harmonics in the breech force could cause im- 
portant dynamic effects. Since the actual breech force is not a simple 
half sine wave, but in fact is a complicated function of time, an investi- 
gation was made to consider the problem further. 

To account for the higher harmonics in the breech-force, the analysis 
decomposes the actual force into a Fourier series with 27 terms. This 
force is applied to the case, and the internal axial force is computed at 
several locations along the case over the duration of the loading. This 
analysis then establishes the contribution of the higher harmonics in the 
breech force to the dynamic stresses in the bar. The analysis is restricted 
to the free-free motor case. 
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ANALYSIS 



In the figure below, AB is the motor case, treated as a one- 
dimensional bar, under the applied breech force P(t). The actual 
displacement of any point of the bar with respect to the X, Y, 
fixed coordinate system is denoted u. This motion consists of two 
parts, a rigid body motion u and a strain-causing deformation u. 
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Then 

u = u + u (1 ) 



If m is the total mass of the bar, then the rigid body 
motion is 

d 2 u = P(t) 



dr 



m 



The differential equation of motion for any point of the bar is 



3 2 0 



c 2 4=0 



( 2 ) 



(3) 



at ax 

where c is the speed of wave propagation, and x is a coordinate 

Z 

axis parallel to x moving with velocity u. The constant c is 
given by 

c = VI 

where E is Young's modulus of elasticity of the bar material and 
p its specific gravity. Combining equations (1) and (3) gives 

^2 (u + u) - c 2 ^ = 0 



(4) 



at 



3x 



2 



3JJ_ 

Substituting for ^2 from (2) yields, 

'in + PLti . c 2 is. = 0 

at 2 " ax 2 

a 2/ \ 2/ , 
With the notation ( ) denoting — , and ( ).. denoting ' 

9x^ tr 3t^ 

the governing differential equation is 



u tt - c u xx 



P(t) 

m 



t > 0, 

0 < x < l 



(5) 



The initial conditions and boundary conditions for the free-free 
system are 

I.C.'s: u(x,o) = o u t (x,o) = o (6) 

B.C.'s: u x (o,t) = u x (£,t) = o (7) 



The solution of equations (5) through (7) gives the deformation of 
any point for any time greater than zero. The stress is then simply 
obtained from Hooke's Law, a = Eu . 

X 

The solution method proceeds by transformation of the given system to 
a system with homogeneous boundary conditions. This is achieved by 
considering a solution in the form 

u(x,t) = v(x,t) + U(x,t). (8) 

Boundary condition equations (7) become 



v x (o,t) = - U x (o,t) 

v x (£,t) = -U x U,t) 

These equations become homogeneous when U(x,t) is taken as 

u(x,t) = (x - ^ ) £ lr 



(9) 



( 10 ) 
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With equations (8) and (10), governing equations (5) through (7) 
become 



v tt - c v xx ' f(x>t) 



( 11 ) 



I.C.'s: v(x,o) = -U(x,o) v t (x,o) = -U t (x,o) (12) 



B.C. 1 s : v (o,t) = o 

A 



v t (£,t) = o 



(13) 



where 



f ( x,t ) " U tt + c2u xx " 

= - (x - d ) Ptt 2P 



04) 



21 ' AE plA 

A Fourier series solution which satisfies boundary conditions (13) is 
assumed for v(x, t) , i.e., 

oo 

~ CM (15) 



V (x , t ) = ^ v n (t) cos 



n=o 



Accordingly, equations (11) and (12) become 



d 2 v (t) 2 

-^2— + “n *„<*> ' f n<») 



(16) 



I.C. 's: v n (o) = o 



dv 



dt 



- (o) = I 



_ 2 | 1 ( x " 21 ) dP(t) 



l 



EA 



dt 



riTrx 

cos -ft- dx 



(17) 



where 



f„(t> = 



2 j ~ / \ riTix . 

f (x, t) cos dx, 



(18) 



and 



rmc 

(l) - — 7 , 

n l 



(19) 
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To this point the analysis is identical to that presented in reference 
1. In that analysis, P(t) was approximated as a half sine wave P 0 sin -nt/a; 
here, the breech force is approximated by the Fourier expansion, 



N 

P(t) = A q + 

k=l 




n 2kri t . n 

I A k cos — + B k 



sin 



2kjTt 1 
a ) 



( 20 ) 



For the linear system, the response due to equation (20) is obtained by 

superposition of individual terms. 

(i) for the uniform force A : 

o 

2P 

From equation (14), f(x,t) = - . Then equation (18) gives 

f = o for all n. The initial conditions are v n = o, and 
v n = o from equation (17), and hence v^(t) = o. Then the 
stress is 

0 ' Eu x - EU x - <' - i> £ < 21) 

This is the stress due to the "rigid body" acceleration for 
constant force. 

(ii) for the cosine terms, A^ cos -^- TT — : 



From equation (14) we obtain. 



2c 

EA£ 



cos 



2k?r t 
a 



Then equations (17) and (18) give 



dt v n<°> = 0 
f n (t > = 

n a n cA 



A^ cos 



2kirt 



( 22 ) 

(23) 
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With equations (22) and (23), the solution of equation (16) 
is 



.. m _ K n a A k f _ 2k7it ^ i 

v „ 1 1 ; — ^ ^ j cos — - — cos a) t y 



n' 



2, 2 n 

» (y„ - 1) 



where 



and 



K = 
n 



Y n = 



21 



? 7 

n a^EA 



2k7T 



Then the stress associated with A, cos - terms is 

K a 

» = Eu x = E(v x + U x ) 



2 

7mA 



(Vn - D 



COS - COS OJ_t 



1 



/ , X , A k 2kxt 

+ < 1 - 7 > r cos — 



sin 



2kirt 



The latter term in equation (27) is the contribution associ 
ated with the "rigid body" acceleration, 
iii) for sine terms, sin 
Equation (14) gives 

x 2 

f ( x , t ) = ( x " 2 I 2ku‘ 



( x " 2 1 ( 

1 EA l 



2c 

EAT 



A. cos 
k a 



From equations (17) and (18), 
d 4k£B k 

k v " (o) 



W ■ ViZ \ sin 

a n LA 



2krr t 



(24) 

(25) 

(26) 

mrx 

l 

(27) 



(28) 

(29) 
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The solution of equation (16) with equations (28) and 
(29) is 



v n (t) 



2K kB.a 
n k 

V 2 



i Y, 



- 2 



- 1 



sin o) n t + 



K a 2 B, 
n k 



ir 2 (Yn " 1 ) 



sin 



2 k it t 

a 



where 

(0 

d = n 

n (n/a ) 



The stress associated with B, sin — - T ~ - is 

k a 




Again the latter term in equation (32) results from the "rigid 
body" acceleration. 

The stress at any point is obtained as the superposition of 
equations (21), (27), and (32). 



(30) 

(31) 

nirx 

l 

(32) 
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RESULTS AND DISCUSSION 



Equations (21), (27), and (32) for the stress have been evaluated 
on the NPS - IBM 360/67 digital computer for specific values of and 
obtained by a Fourier decomposition of the given breech force. This 
decomposition is done by the IBM subroutine FORIT on the computer. The 
present analysis resolved the breech force into 26 sine and cosine har- 
monics, i.e., 

26 

P(t)=A o + I f (33) 

k=l 1 

The magnitude of the harmonic coefficients A^ and B^ depend upon the 
nature of the breech force. For this analysis, the force given in Fig. 

1 of Ref. 2 was used. Approximation of that force by equation (33) leads 
to 
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-.0940 
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-.0552 
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-.8107 




17 
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-.1162 


.1077 




18 


-.2053 
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5 


.1089 


-.1257 
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-.1821 


.0827 
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-.2036 


.3394 
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-.2162 


.0815 
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-.3478 


-.0886 




21 


-.2723 


.0286 
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-.3001 


-.1664 




22 


-.2164 


-.0294 
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-.2229 


-.3071 




23 


-.1599 


.0250 


10 


.1728 


-.2293 




24 


-.2146 


.0638 


11 


.2260 


.1067 




25 


-.2589 


-.0065 


12 


-.0425 


.3243 




26 


-.2162 


.0034 


13 


-.3157 


.2287 
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A plot of equation (33) with the above values of A^ and is 
given in Fig. 1 for a .026 second time duration of breech force. 

Figures 2 through 5 are the results for the stress in the bar due to 
the breech force shown in Fig. 1 as a function of time for j ratios 
of 0, 1/4, 1/2, and 3/4, respectively. (The analysis is for the 5 in. 
steel shell of 0.19 in. thickness). Comparing Fig. i with Figs. 2 
through 5 reveals that the internal force appears as a duplicate of 
the breech force, with the stress curves scaled versions of the force 
input. In Figs. 6 through 8, the contributing harmonics are individu- 
ally displayed. The elastic dynamic effects are the small oscillations 
about each harmonic itself. 

Figures 10 through are the results of the analysis for a 
time duration of 0.0052 seconds for the brech force. This is a re- 
duction in the actual duration by a factor of five. This sample prob- 
lem was considered in order to assess the influence of the magnitude 
of the force duration. The elastic dynamic effects can be clearly 
seen in these figures. Note that the magnitude of stress has been 
significantly changed by the dynamic effects. 

CONCLUSIONS 

The results of this analysis (Figs. 1-5) show that dynamic 
effects are negligible for the breech force of 0.026 seconds time 
duration of Figure 1, Reference 2. Hence, when the propellant is 
not bonded to the motor case, the case responds to the breech force 
as if it were applied in a static sense, i.e., the maximum internal 
axial load in the case can be computed on the basis of the mass 
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distribution and the maximum rigid body acceleration of the rocket. 
This is the same conclusion arrived at in Ref. 1 based on a simple 
one-half sine wave pulse; hence, accounting for the more complicated 
actual breech pressure pulse confirmed the original conclusion. 

Further investigation (Figs. 10-14) verified the dependence 
of the dynamic effects on the presence of significant higher harmonic 
terms in the breech force. Such harmonics may be the result of shor- 
tening the time duration of loading, or by the presence of a local 
rapid perturbation of force on an otherwise acceptable smooth breech 
force input. 
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